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' 1- Introduction 



Many years ago, in 1967[1], H. R. Lewis has shown that there is a conserved quantity, 
that will be indicated by 1, associated with the time dependent harmonic oscillator (TDHO) 
with frequency uj{t), given by: 

/ = lliqa - aqf + {^f] , (1.1) 

where q and a obey, respectively the equations: 

q + uj\t) g = , a + Lo\t) a = 4 . (1.2,3) 
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On the other hand, as the above expressions have also been obtained by V. P. Er- 
makov [2] in 1880, the invariants determination of time dependent physical systems is also 
known as the Ermakov-Lewis problem. So, considerable efforts have been devoted to 
solve this problem and its generalizations, in the last thirty years, and in many works have 
been published on these subjects [3,4]. 

In the present work we investigate the existence of these invariants for the one- 
dimensional non-linear Gross-Pitaevskii equation with the potencial V{x, t) given by: 

V(x, t) = Im u^it) , (1.4) 

which is the time dependent harmonic oscillator potential. 
2. Gross-Pitaeviskii Equation 

Em 1961 [5, 6], E. P. Gross and, independently, L. P. Pitaevskii proposed a non-hnear 
Schrodinger equation to represent time dependent physical systems, given by: 

. ^ ^ - ^ + i m u^\t) ^(x, t) + g\ i^{x, t) f , (2.1) 

where ip{x, t) is a wavefunction and ^ is a constant, constante. 

Writing the wavefuncition ip{x, t) in the polar form, defined by the Madelung-Bohm 
transformation[7,8], we get: 

il){x, t) = 0(a;, t) e' *) , (2.2) 

where S{x , t) is the classical action and (p{x, t) will be defined in what follows. 

Substituting Eq.(2.2) into Eq.(2.1) and taking the real and imaginary parts of the 
resulting equation, we get [9]: 

t + = ' (2-3) 

^ + ^ + u;^{t) x = - ^ £ (K. + Vgp) , (2.4) 

where: 

p{x, t) = (j)'^{x, t) , (2.5) (quantum mass density) 
Vgu(x, t) = J , (2.6) (quantum velocity) 

Vqu{x, = - ^ 7^ , (2.7) (Bohm quantum potential) 



and 
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Vgp = ^ P ■ (2-8) (Gross-Pitaevskii potential) 

In order to integrate Eq.(2.4) let us assume that the expected value of quantum 
force is equal to zero for all times t, that is, 

<^> ^0 ^ ^U = lit) > < ^ > = q{t) . (2.9a-c) 
In this way, we can write Eq.(2.4) into two parts: 

^ + v,^^ + co'x ^ k{t) [x - q{t)] , (2.10) 

- ^ I: + VGP) = £{S^j=/-^-ip) = m [x - q{t)] . (2.11) 

Performing the differentiations indicated in Eq.(2.11) we get, 

[pi^~^lfi^ + ^ + m If ^ ^(^) ~ 5(^)1 ■ (2-12) 

To integrate Eq.(2.12) it is necessary to known the initial condition for p{x, t). Let 
us assume that for t = the physical system is represented by a normalized Gaussian wave 
packet, centered at q{0), that is, 

p{x, 0) = [tt a(0)]- V2 e ^(5^ = e" - , (2.13) 
A = TT (7(0) , S = a; - g(0) , C = (7(0) . (2.14-16) 
Since Eq.(2.13) is a particular solution of Eq.(2.12), we must have: 

n"^ r 1 d^pjx, 0) _ 2 dpjx, 0) 0) , 1 ( dpjx, 0) n31 , fl. QpCa:, 0) ^ 

4 m2 \-p{x, 0) 9x3 p2(j.^ 0) dx dx^ ' p{x, 0)^ ^ dx ^ i ' m dx 

= k{t) [x - q{t)] . (2.17) 
Performing the differentiation indicated above, Eq.(2.17) becomes: 

] [x - qm = m) [x - q(Q)] ^ 

A;(0) = ^ ] ■ (2.18) 

Comparing Eq.(2.18) with the Eqs.(2.12) and (2.13), by analogy we get. 



kit) = ^ ] , (2.19) 

t) = [tt V2 g- "-(?)*^'' , (2.20) 

where, 

^ ' ait) m vVa(t) 

Taking into account Eqs.(2.19-21), let us perform the foUowing differentiations, re- 
membering that t and x are independent variables: 

t = - H P + P ^ + "-^Q] , (2-22) 

If = - P , (2.23) 

^From Eq.(2.3) and Eqs. (2.22-23), results 

^ - = ^ - ^ - ,)2 _ liE^ q . (2.24) 

Defining 

p{x, t) ^ - 2 ^ , (2.25) 

r{x, t) = f-^ - ^{x - qf - q , (2.26) 

Eq.(2.24) becomes, 

^ + p{x, t) = r{x, t) , (2.27) 
which can be integrated, giving: 

v,u = liJrudx + c{t)] , (2.28) 

where, 

u = exp (/ p dx) . (2.29) 
Using Eqs.(2.20,28,29), the function u given by Eq.(2.29) is written us[9]: 
u ^ exp[![- ^^2^] ax ) = (tt (7)V2 p . (2.30) 
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In this way, defining / = r u x, and using Eqs. (2.26,30) we obtain: 

/ = f r u dx — 

- I[f-^ - ^{x - qy - '-^q]{nay/'pdx - 

I ^ h - h, (2.31) 

where, 

h = I[f^ - ^{x - q)'] (tt a)V2 p dx , (2.32) 

and 

h = li^^^q] {nay/'pdx. (2.33) 

To integrate Eq.(2.32) it is necessary, first to perform the differentiation [9] shown 
bellow, where Eq.(2.23) is used: 

£[f-A^- <I)P] - if-. - ^^]P- (2-34) 
Inserting Eq.(2.34) into Eq.(2.32), results 

h = (tt ay/' p (^) {x - q). (2.35) 

Similarly, to calculate /2, seen in Eq.(2.33), we need to use Eq.(2.23) obtaining, 

/2 = - (tt (7)1/2 ^ p _ (2.36) 

Substituting Eqs. (2.35-36) into Eq.(2.31) we see that, teremos: 

/ = (tt a)V2 p [_£_ _ 5) + ^] . (2.37) 

Remembering that the quantum velocity Vq^ is defined by Eq.(2.28) and using Eqs. (2.30,37) 
we verify that Vqu can be written as: 

v,n = f^{x - q) + q + ^^^^ . (2.38) 

Assuming that the mass density p — > when | x | — > oo we verify that the parameter 
c{t) must be equal to zero. Consequently, Vqu becomes, 

^) = 2^ - 9(0] + m ■ (2-39) 
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Using the above Eq.(2.39) we calculate the following differentiations, remembering 
that t and x as independent variables: 

- iif-A^ - <i) + ^] - 

= A - ^) - ^ - ^) - A ^ + ^ ■ (2-40) 

dVnu 



dx dx 



(x - + g] = ^ . (2.41) 



Now, adding the factor a;^ q to Eqs. (2.39-41), we see that Eq.(2.10) becomes, 

2^ - - 1^ - - 2^ 9 + ? + [l^ - + 9] 2^ + 

+ ^2 ^ ^ ^2 g _ ^2 ^ ^ (_|L, _ _^) _ ^) ^ 

- + a;2 + ^3 ] _ ^) + ^- + ^2 . (2.42) 

To satisfy Eq.(2.42), the following conditions must be obeyed: 

^ _ + + = , (2.43) 

g + g = . (2.44) 

Putting 

a = A «2 ^ (2.45) 

we obtain, 

a ^ ^ aa , ^ = ^ [(a)^ + a a] . (2.46-47) 
Inserting Eqs.(2.45-47) into Eq.(2.43): 

2 n [(d)^ + « g] _ " I , ,,2 _ n ^ 



a + a;2 ^ ^ j_ _ (2.48) 
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Note that, although the above Eq.(2.48) is formally identical to that obtained by 
Ermakov [see Eq.(1.3)], they are different due to the Planck's constant h which appears in 
Eq.(2.48). This is the same kind of difference found between the classical wave equation 
(d'Alembertian) and the quantum Schrodinger wave equation. 

Finally, eliminating the factor o;^ into Eqs.(2.44) and (2.48) we get, 

y m V m 



_2L _ 2 .'V q 



I (d g - qa) 

{a q - q a) {a q - q a) = - ^ J'^^ ) {a q - q a) 



|[i {aq - qaf] = - ^|m + ^^^-^ 



|[i {aq - qaf + \ {^f 



2 .9 q A [a.] 



where. 



§ = -I (1) , (2.49) 

at / h at ^a' ' ^ ' 



I ^ \ \{aq - qaf (f )^] , (2.50) 



which represents the Ermakov-Lewis-Schrodinger invariant of the time dependent har- 
monic oscillator (TDifO)[9]. In conclusion, we have shown that the Gross-Pitaesvskii 
equation has not an Ermakov-Lewis invariant for the TDHO. 
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